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Unity  transmittance  at  an  interface  between  bulk  media  is  quite  common  for  polarized  electromagnetic 
waves  incident  at  the  Brewster  angle,  but  it  is  rarely  observed  for  sound  waves  at  any  angle  of  incidence.  In 
the  following,  we  theoretically  and  experimentally  demonstrate  an  acoustic  metamaterial  possessing  a 
Brewster-like  angle  that  is  completely  transparent  to  sound  waves  over  an  ultra-broadband  frequency  range 
with  >100%  bandwidth.  The  metamaterial,  consisting  of  a  hard  metal  with  sub  wavelength  apertures, 
provides  a  surface  impedance  matching  mechanism  that  can  be  arbitrarily  tailored  to  specific  media.  The 
nonresonant  nature  of  the  impedance  matching  effectively  decouples  the  front  and  back  surfaces  of  the 
metamaterial  allowing  one  to  independently  tailor  the  acoustic  impedance  at  each  interface.  On  the 
contrary,  traditional  methods  for  acoustic  impedance  matching,  for  example  in  medical  imaging,  rely  on 
resonant  tunneling  through  a  thin  antireflection  layer,  which  is  inherently  narrowband  and  angle  specific. 

The  concept  of  metamaterials  was  originally  introduced  for  electromagnetic  waves  and  received  significant 
attention  over  the  last  decade,  due  to  unconventional  interaction  of  these  man-made  materials  with  waves,  as 
in  the  case  of  negative  refraction  materials \  More  recently,  the  metamaterial  concept  has  been  extended  to 
acoustic  waves  in  a  variety  of  scenarios  of  interest  such  as  acoustic  clocking^,  superlensing^  and  sound  focusing 
and  confinement^.  Several  challenges  still  hold  in  applying  metamaterials  to  realistic  devices,  in  particular  when 
low  losses  and  wide  bandwidths  of  operation  are  desired,  as  most  of  their  exotic  features  are  based  on  resonant 
inclusions. 

One  particular  class  of  metamaterials  that  has  received  special  attention  in  recent  years  is  an  array  of  sub- 
wavelength  apertures  in  an  opaque  screen,  producing  spectral  bands  of  extraordinary  transmission.  In  other 
words,  more  energy  passes  trough  these  sub  wavelength  apertures  than  one  would  expect  from  simple  geometrical 
considerations  based  on  the  aperture  to  period  ratio.  At  first,  this  phenomenon  has  been  predicted  and  experi¬ 
mentally  observed  in  the  electromagnetic  domain,  where  it  is  known  as  extraordinary  optical  transmission 
(EOT)^“^.  Lately,  it  has  been  extended  to  acoustic  waves,  as  extraordinary  acoustic  transmission  (EAT)^““. 
Usually  these  phenomena  are  achieved  under  some  form  of  resonant  condition,  and  are  therefore  limited  in 
bandwidth  and  lossy,  particularly  in  the  case  of  very  subwavelength  apertures.  In  this  paper,  we  demonstrate  a 
different  way  of  squeezing  and  tunneling  acoustic  waves  in  ultranarrow  apertures,  overcoming  several  of  the 
conventional  EAT  limitations.  To  this  goal,  we  show  a  way  to  manipulate  the  effective  constitutive  properties 
(density  peff  and  sound  velocity  Cgff)  of  an  acoustic  metamaterial  mostly  composed  of  impenetrable  hard  materi¬ 
als,  in  order  to  realize  broadband  impedance  matching.  As  an  example  of  the  powerful  versatility  of  this  imped¬ 
ance  matching  technique,  in  the  following  we  demonstrate  acoustic  matching  of  the  metamaterial  with  air.  Noise 
in  the  everyday  environment  is  ubiquitous,  primarily  due  to  the  large  impedance  mismatch  between  air  and  solid 
materials.  Most  rigid  acoustic  materials  (glass,  steel,  etc.)  have  a  large  acoustic  impedance  compared  to  air,  which 
makes  them  practically  impenetrable  to  sound  waves.  This  is  related  to  the  fact  that  the  angle  of  acoustic 
intromission^^,  i.e.,  the  angle  at  which  the  reflection  is  zero  for  an  interface  between  air  (medium  1)  and  a  generic 
material  (medium  2)  with  effective  density  and  effective  sound  velocity  c^ff ,  is  rarely  found  in  conventional 
materials.  When  it  exists,  the  intromission  angle  {Sj)  must  satisfy  the  following  expressions^ 


sin  Si  = 

\ 


(1) 
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where  po  and  Cq  are  respectively  the  density  and  sound  velocity 
in  medium  1  (air  in  our  case).  The  angle  for  unity  transmittance 
is  similar  to  the  Brewster  angle  for  electromagnetic  waves  and, 
although  the  underlying  mechanisms  for  perfect  transmittance  of 
acoustic  and  electromagnetic  waves  are  very  different,  they  can  both 
be  understood  as  an  impedance  matching  phenomenon. 

It  is  seen  in  Eq.  (1)  that  the  intromission  condition  may  be  met 
when  either  PeffJ Po>Cq/c I  or  Pef^Po^ColCef<l,  which  cannot  be 
satisfied  by  most  ordinary  bulk  materials,  due  to  the  fact  that  the 
speed  of  sound  generally  increases  with  density.  One  instance  of  a 
naturally  occurring  intromission  angle  is  in  ocean  acoustics  at  the 
interface  of  water  and  the  sea  floor^^  For  a  muddy  bottom  composed 
of  water  saturated  silt,  the  speed  of  sound  is  lower  than  water  but  the 
density  is  higher  than  water.  This  is  the  intromission  condition  that 
we  describe  below;  a  metamaterial  having  a  density  higher  than  air 
but  a  speed  of  sound  lower  than  air. 

In  the  following,  we  theoretically  predict  and  experimentally  dem¬ 
onstrate  that  an  effective  angle  of  intromission  may  be  introduced  at 
will  in  a  metamaterial  formed  by  an  array  of  sub -wavelength  aper¬ 
tures  in  a  thick,  acoustically  hard  screen.  At  this  angle,  nearly  total 
transmission  spanning  an  ultrabroad  range  of  acoustic  frequencies  is 
obtained.  Most  EAT  phenomena  rely  on  resonance  effects  that  are 
inherently  narrow-band,  and  for  which  large  transmission  is  usually 
hindered  by  absorption  losses  as  the  aperture  size  becomes  much 
smaller  than  the  wavelength.  In  addition,  resonant  EAT  frequencies 
strongly  depend  on  the  screen  thickness  and  the  array  period,  relying 
on  the  resonant  coupling  between  entrance  and  exit  of  each  aperture. 
Here,  on  the  other  hand,  the  nonresonant  nature  of  the  intromission 
tunneling  separately  arises  at  the  entrance  and  exit  of  the  screen, 
implying  that  is  very  robust  to  absorption  or  energy  extraction  for 
harvesting  purposes  within  the  apertures,  and  may  be  tailored  to 
impedance  match  dissimilar  media.  The  impedance  matching  is 
essentially  frequency  independent  and  may  be  tailored  by  the  geo¬ 
metry  of  the  metamaterial,  as  we  describe  in  the  following. 


Results 

In  order  to  provide  some  physical  insights  into  this  phenomenon, 
we  refer  to  the  geometries  described  in  Fig.  1:  an  acoustic,  time- 
harmonic  c  plane  wave  with  wave-number  ko=(JL>/co  is  incident 
at  an  arbitrary  angle  {S,(p)  on  an  acoustically  hard  screen  of  thickness 
/,  corrugated  by  slits  of  width  w  and  period  d;  we  also  consider  the 
case  of  2-D  apertures  with  periods  and  dy ,  area  of  the  elementary 
cell  Si  =  dxdy  and  area  of  the  apertures  S2  . 


Figure  1  |  Sketch  of  the  considered  geometries.  On  the  left  a  1-D  acoustic 
grating  of  period  d,  slit  aperture  w  and  thickness  /;  on  the  right  a  2-D 
acoustic  grating  of  periods  d^  and  dy ,  area  of  the  elementary  cell  Sj  =  d^dy 
and  area  of  the  aperture  S2.  Both  structures  are  excited  by  an  acoustic 
plane- wave  incident  at  an  angle  3  with  respect  to  the  z-axis  and  an  angles 
in  the  x-y  plane. 


Here  we  suppose  that  the  grating  periodicities  are  smaller  than  the 
wavelength  of  the  incident  sound  wave  (X>2dx,  X>2dy)  so  that  all 
the  diffraction  orders  except  the  zero-th  are  evanescent.  Under  this 
condition  reflected  and  transmitted  pressure  fields  may  be  respect- 

ively  written  as  =  +  for 

z<0  and  =  for  z>Z  where  r  and  t 

are  the  reflection  and  transmission  coefficients,  respectively,  and 
ky  are  the  transverse  components  of  the  wavevector,  which  may 
be  written  in  terms  of  the  angles  3  and  (p  (see  Fig.  1)  as: 
kx=kosm{S)cos{(p)  and  k^=A:osin(,9)sin(^).  Inside  each  aperture, 
only  the  (0,0)  acoustic  mode  can  propagate,  with  general  expression 
Pwg(z)  =  -\-Be~ .  Unlike  electromagnetic  waveguides, 

there  is  no  cutoff  wavelength  for  a  2-D  acoustic  aperture.  In  addition, 
due  to  momentum  conservation,  each  channel  is  separately  excited 
with  a  phase  that  matches  the  impinging  excitation  at  the  location 
(xo,yo)  of  aperture,  implying  that  we  may  write  for  the  pressure 
fields  =  +  for  0<z<l,  where  A 

and  B  are  the  unknown  amplitudes  of  the  forward  and  backward 
propagating  modes.  By  imposing  the  continuity  of  the  pressure  field 
P  and  the  normal  component  of  the  volume  velocity  S(ipo(jj)~^dPldz 
at  z=0  and  z=l,  the  transmission  takes  the  following  expression: 


t  = 


4 


^acl 


,  (2) 


where,  according  to  the  notation  in  Ref  we  have  defined  the  acous¬ 
tic  impedance  of  the  impinging  wave  on  the  surface  Si  as  Zaci=PoCol 
(SiCOSi9),  and  the  acoustic  impedance  of  the  guided  mode  inside  the 
channels  of  cross  section  S2  as  Zac2  =  PoColS2.  For  2-D  apertures  (Fig.  1 
right  panel)  Si=dxdy  is  the  area  of  the  elementary  cell  and  S2  is  the 
area  of  the  aperture,  while  for  1-D  slits  (Fig.  1  left  panel )  Si  =  lm'd 
and  S2=lm-w  are  the  areas  normalized  to  i  m  along  the  y-axis. 
Resonant  EAT  (|t|^  =  l)  is  obtained  at  the  FP  resonances  kol=n7i, 
with  n  being  an  integer,  consistent  with  the  findings  in  Refs.^““.  Yet, 
Eq.  (2)  admits  another  peculiar  condition  for  unitary  transmission, 
which  arises  when  Z^^i—^acz-For  normal  incidence  this  condition 
cannot  be  met,  due  to  the  difference  in  cross  sectional  areas,  but,  by 
increasing  the  angle  of  incidencei9,  can  be  smoothly  increased 
until  it  matches  the  value  of  Z^^z  at  the  angle  Sf. 

cos  81  =  82/ Si,  (3) 


which  represents  the  intromission  angle  for  an  acoustic  metamaterial. 
Although  rare  for  ordinary  bulk  materials an  intromission  angle 
can  always  be  obtained  for  the  metamaterial  shown  in  Fig.  1  and  the 
angle  depends  only  on  the  geometry,  i.e.  the  percentage  of  the  open 
area  of  the  screen.  This  anomalous  matching  phenomenon  is  totally 
independent  on  the  grating  thickness,  since  the  effective  channel 
impedance  matches  the  impinging  wave  at  each  interface,  producing 
anomalous  total  tunneling  and  energy  squeezing  through  each  chan¬ 
nel.  This  strikingly  simple  formula  captures  to  a  large  degree  a  novel 
tunneling  phenomenon  which  may  be  regarded  as  the  acoustic  ana¬ 
log  of  the  “plasmonic  Brewster  angle”  introduced  in  Ref  for  1-D 
metallic  gratings  excited  by  electromagnetic  waves.  Due  to  the 
absence  of  modal  cut-off  in  2-D  acoustic  waveguides,  this  condition 
interestingly  holds  also  for  2-D  gratings  of  arbitrary  cross-section  in 
the  case  of  acoustic  waves.  Eq.  (2)  is  consistent  with  a  rigorous 
homogenization  of  such  acoustic  grating  as  a  homogeneous  acoustic 
metamaterial  with  effective  spatially-dispersive  parameters: 


sm^{9),p,ff  =  PoSi/S2. 


(4) 


Indeed,  the  acoustic  intromission  angle  as  defined  in  Eq.  (1)  is  given 
by  Eq.  (3).  We  notice  that  in  our  model  there  is  no  dependence  on  the 
azimuthal  angle  (p,  as  one  may  expect  given  the  scalar  nature  of  the 
pressure  fields.  The  advantage  of  the  nonresonant  mechanism  is 
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Figure  2  |  Experimental  set-up.  Schematic  drawing  of  the  experimental 
set-up  and  images  of  the  1-D  and  2-D  metamaterials. 


that  acoustic  energy  collectors  and  sensors  can  be  placed  inside  the 
apertures  without  harming  the  impedance  matching  and  funneling 
at  the  surface.  What  is  going  on  inside  the  apertures  doesn’t  affect  the 
collection  efficiency  as  it  would  for  a  resonant  mechanism!  In  addi¬ 
tion,  the  waveguides  can  reroute  the  acoustic  energy  into  a  single 
waveguide  for  energy  collection  onto  one  device  instead  of  many.  By 
tapering  the  cross  section  area  of  the  waveguides  with  depth  into  the 


metamaterial  it  is  possible  to  change  the  open  area  of  the  screen  for 
the  incident  and  exit  surfaces  independently,  thereby  achieving 
impedance  matching  to  two  acoustic  media  with  different  values  of 
impedance. 

In  order  to  verify  the  basic  prediction  of  acoustic  impedance 
matching  to  air,  we  have  fabricated  1-D  and  2-D  subwavelength 
aperture  arrays  in  a  thick  aluminum  plate  and  measured  the  trans¬ 
mission  as  a  function  of  incidence  angle  (0)  and  acoustic  frequency. 

Acoustic  measurements  were  performed  in  an  acoustic  anechoic 
chamber  in  the  frequency  range  4-20  kHz:  a  sketch  of  the  experi¬ 
mental  set-up  and  pictures  of  the  samples  are  shown  in  Figure  2.  The 
sample  was  placed  on  a  rotational  stage  and  the  sound  pressure 
transmission  measurements  were  performed  from  normal  incidence 
to  nearly  grazing  angle.  In  Fig.  3  we  show  the  experimental  results  for 
the  case  of  a  1-D  structure  and  compare  them  with  our  theoretical 
homogenized  model  (2-4). 

The  upper  panel  of  Fig.  3  corresponds  to  the  1-D  aluminum  sam¬ 
ple  with  dimensions  1=25.4  mniy  d=4.275  mm  and  w=  1.1  mm.  Our 
analytical  model  captures  with  remarkable  accuracy  the  fundamental 
physical  mechanisms  behind  the  transmission  resonances  of  the 
grating:  at  normal  incidence  (^9=0),  typical  EAT  peaks  based  on 
FP  resonances  are  visible.  These  resonances  are  the  horizontal  bands 
in  the  plots,  since  this  EAT  mechanism  is  inherently  narrow-band, 
but  weakly  dependent  on  the  incidence  angle.  In  contrast,  the  intro¬ 
mission  angle  transmission  arises  as  a  vertical  band,  confirming  weak 
dependence  on  frequency,  but  selectivity  to  the  transmission  angle. 
One  may  wonder  how  narrow  can  the  slits  be  made  and  still  realize 
high,  ultra-broadband  transmission  at  the  intromission  angle,  con- 
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Figure  3  |  Transmittance  for  the  1-D  samples.  Angular  power  transmission  spectra  for  two  different  1-D  gratings:  (upper  panel)  /=25.4  mm, 
d=  4.275  mm  and  w=  1.1  mm;  (lower  panel)  /=25.4  mm  d=  3.575  mm  and  w=400  pm.  Transmission  bands  independent  of  angle  are  due  to  resonant 
Fabry-Perot  modes  and  transmission  bands  independent  of  frequency  are  due  to  the  non-resonant  intromission  process.  Open  area  of  the  screen  for  the 
two  cases  are  26%  (upper)  and  11%  (lower). 
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Figure  4  |  Transmittance  for  the  2-D  sample  and  localization  of  the  pressure  field  inside  the  channels.  Angular  power  transmission  spectra  for  a  2-D 
hole  array  aluminum  grating  with  1=25.4  mm,  dx=  dy=  7.94  mm  and  hole  radius  r=2  mm.  Pressure  wave  numerically  calculated  for  the  second  FP 
resonance  at  v=  12.75  KHz  in  the  case  of  normal  incidence  and  at  the  intromission  angle  as  indicated  in  the  figures.  The  open  area  of  the  screen  is  20%. 


sidering  that  conventional  EAT  phenomena  are  very  sensitive  to 
presence  of  losses  in  the  channels,  in  particular  for  large  field  con¬ 
finements  in  small  apertures  and  slits.  In  our  homogenized  model 
(2-4),  we  have  neglected  the  finite  frictional  losses  and  field  penetra¬ 
tion  in  the  metal,  although  they  are  inevitable  in  our  experiments.  In 
order  to  assess  the  robustness  of  the  broadband  intromission  tunnel¬ 
ing  phenomenon  against  frictional  losses,  we  have  fabricated 
and  tested  a  second  1-D  sample  (see  lower  panel  of  Fig.  3)  with 
same  thickness  1=25.4  mm,  but  different  period  and  slit  width 
d= 3.575  mm,  w=400  /im.The  open  area  of  the  screen  is  11%  . 

Our  analytical  model  has  been  extended  to  account  for  small 
frictional  losses  by  introducing  a  perturbative  dissipation  factor  in 
the  wave-vector,  which  becomes^®  ko[l  -\-  (2w))],  where  is 

the  thickness  of  a  viscous  boundary  layer.  For  air  at  20°C,  the  ex¬ 
pression  for  the  viscous  boundary  layer^^  is:  d^^^  =  (2.19/Vv)  mm 
where  v  is  the  frequency  in  Hz  and  is  —20  pm  thick  at  these  fre¬ 
quencies.  At  normal  incidence,  the  FP  transmission  bands  are 
severely  affected  by  the  presence  of  losses  in  this  second  scenario, 
both  in  our  analytical  results  and  in  our  experiment,  when  compared 
to  the  case  with  1.1  mm  slits.  The  intromission  tunneling  is  however 
significantly  more  robust,  providing  one  order  of  magnitude  larger 
transmission  than  at  the  conventional  FP  resonances.  This  robust¬ 
ness  is  a  symptom  of  the  non-resonant  nature  of  the  intromission 
tunneling.  In  conventional  FP  tunneling,  the  pressure  field  under¬ 
goes  a  series  of  multiple  reflections  at  the  highly  mismatched 
entrance  and  exit  faces  of  the  grating.  This  effectively  forms  a  res¬ 
onant  cavity  with  large  stored  energy,  thereby  enhancing  the  fric¬ 
tional  losses.  In  contrast,  broadband  intromission  tunneling  and 
energy  squeezing  is  achieved  due  to  impedance  matching,  for  which 
the  impinging  wave,  independent  of  the  frequency  of  excitation, 
‘feels’  the  same  impedance  in  vacuum  and  inside  the  slits,  avoiding 
any  reflection  at  the  entrance  and  exit  faces  of  the  grating,  regardless 


of  the  grating  thickness  1.  As  the  filling  ratio  gets  smaller,  the  sound 
squeezing  increases  and  the  corresponding  angular  spread  of  the 
intromission  transmission  band  gets  narrower,  moving  towards  stee¬ 
per  angles.  A  clear  vertical  high  transmission  band  at  the  intromis¬ 
sion  angle  is  seen  spanning  the  entire  frequency  range,  even  in  this 
extreme  scenario,  highlighting  the  robustness  of  this  concept  in 
terms  of  its  practical  applications. 

We  now  analyze  the  2-D  sample  which  may  be  considered  even 
more  interesting  for  energy  concentration,  harvesting  and  focusing 
purposes.  In  Fig.  4  we  compare  theory  and  experiments  for  2-D 
apertures.  Even  in  this  2-D  geometry,  a  clear  broadband  intromission 
band  is  observed  in  our  experiment  (right  figure  of  the  upper  panel), 
well  captured  by  our  analytical  model  (left  figure  of  the  upper  panel), 
confirming  that  the  grating  still  retains  the  main  physical  properties 
of  a  metamaterial  (4^  We  notice  a  slight  shift  of  the  measured  spec¬ 
tral  position  of  the  FP  transmission  bands  compared  to  our  analytical 
model;  a  corresponding  small  shift  is  also  observed  in  the  angular 
position  of  the  intromission  band.  Both  effects  are  attributed  to 
having  neglected  in  our  analytical  model  the  reactive  fields  associated 
with  the  evanescent  higher  diffraction  orders  excited  at  the  grating 
interfaces,  which  may  be  taken  into  account  by  considering  a  slightly 
longer  channel  length  and  a  slight  difference  in  the  input  imped¬ 
ance^^.  Due  to  the  large  abruption  in  this  2-D  scenario,  these  effects 
are  expected  to  play  a  somewhat  more  relevant  role  than  in  the  1-D 
case.  In  order  to  further  highlight  the  anomalous  features  of  this 
tunneling  phenomenon,  in  the  lower  panel  we  show  the  numerical 
calculation^^  of  the  pressure  wave  outside  and  inside  the  channel  at 
the  frequency  v  =  12.75  KHz,  which  represents  the  frequency  of  the 
second  FP  resonance  of  the  grating.  We  compare  the  case  of  normal 
incidence  excitation  to  the  case  of  excitation  at  the  intromission 
angle.  It  is  evident  how,  despite  both  excitations  ensure  very  high 
transmission  levels,  the  field  distributions  inside  the  channels  are 
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very  different.  For  normal  incidence,  total  transmission  is  obtained 
by  exploiting  a  strong  FP  resonance  established  via  the  large  mis¬ 
match  at  the  entrance  and  exit  faces,  which  forms  a  typical  standing 
wave  distribution  inside  the  channel.  In  contrast,  at  the  intromission 
angle  tunneling  is  based  on  impedance  matching:  the  wave  does  not 
experience  any  change  in  effective  impedance  and  it  can  tunnel 
through  the  channel  with  nearly  uniform  amplitude,  no  amplifica¬ 
tion  and  minimum  reflections  and  standing- wave  ratio.  The  normal 
component  of  the  velocity  (not  shown  here)  is  the  one  being  amp¬ 
lified  in  the  channel,  due  to  the  continuity  of  the  total  volume 
velocity,  but  similar  differences  are  seen  between  incidence  in  the 
normal  direction  and  at  the  intromission  angle  in  terms  of  uniform 
distribution  and  absence  of  standing  waves  in  the  channel.  We  also 
note  that  as  long  as  the  field  is  tuned  on  a  FP  resonance,  by  increasing 
the  incident  angle  it  gradually  passes  from  a  resonant  channeling 
regime  to  a  non-resonant  channeling  regime  with  no  drop  in  the 
overall  transmittance. 

Discussion 

These  findings  may  have  important  applications  for  sub -wavelength 
acoustic  imaging,  acoustic  field  concentration,  harvesting,  and 
impedance  matching  for  multiple  applications  including  medical 
diagnostics^^  and  noise  controP^.  For  example,  2-D  acoustic  meta¬ 
materials  have  been  explored  for  deep -sub wavelength  imaging^^  by 
exploiting  their  FP  resonances,  also  an  acoustic  broad-band  hyper¬ 
lens  has  been  demonstrated^^.  Our  results  suggest  that  similar  appli¬ 
cations  could  be  devised  at  the  intromission  angle  of  the  grating,  with 
exceptional  broadband  characteristics.  We  have  also  found  that,  even 
by  reducing  the  slit  apertures,  the  intromission  band  remains  much 
less  affected  by  the  frictional  losses  with  respect  to  the  FP  resonances, 
and  thus  gives  the  possibility  to  funnel  sound  through  extreme  sub- 
wavelength  apertures.  This  effect  is,  again,  a  symptom  of  the  different 
physical  mechanisms  involved  in  classic  EAT  based  on  FP- 
resonances  compared  to  the  nonresonant  mechanism  at  the  intro¬ 
mission  angle  that  we  have  elucidated  in  this  work.  Even  in  the 
extreme  case  in  which  the  acoustic  energy  funneled  into  the  2-D 
apertures  would  be  harvested  or  converted  into  other  energy  forms, 
for  which  the  absence  of  a  standing- wave  resonance  would  suppress 
any  form  of  FP  tunneling,  the  matching  phenomenon  described  here 
would  be  preserved.  This  important  property  suggests,  for  example, 
the  application  of  these  concepts  to  realize  broadband  acoustic  to 
electric  energy  conversion  through  piezoelectric  materials  using 
micromechanical  systems  (MEMS)  technology.  Finally,  the  results 
presented  here  for  impedance  matching  in  air,  indicate  future  pos¬ 
sibilities  for  liquid  and  solid  media. 

Methods 

1.  Sample  preparation.  The  ID  metallic  sample  was  constructed  using  ~300 
extruded  aluminum  flat  bars  0.3175  cm  by  2.54  cm.  The  sample  size  was  ~0.5  m 
high,  ~  1  m  wide  and  2.54  cm  thick,  much  larger  than  the  incident  beam  size,  in  order 
to  minimize  any  diffraction  effects  around  the  sample,  particularly  for  large  incidence 
angles.  The  grating  was  held  together  by  metal  rods  through  the  bars  at  the  top  and 
bottom  of  the  sample  and  the  aperture  width  was  fixed  by  inserting  shim  spacers  near 
the  rods  of  1 . 1  mm  and  400  microns.  The  2D  hole  array  structure  was  built  by  drilling 
4  mm  holes  into  a  2.54  cm  thick  aluminum  plate  with  dimensions  of  9 1  cm  by  46  cm. 

2.  Angular  acoustic  transmission  measurements.  The  acoustic  angular 
transmission  measurements  were  performed  by  building  a  small  anechoic  chamber 
using  sound  absorbing  foam.  The  sample  was  placed  on  a  rotational  stage  in  the 
middle  of  the  chamber  and  the  distance  between  the  driver  and  microphone  was 
190  cm.  A  Bruel  &  Kjaer  pulse  system  was  used  together  with  a  50  kHz  Crwon  xTi 
1000  amplifier,  Bruel  &  Kjaer  4961  free  field  microphone  and  BMS  4512ND  planar 
wave  drive  of  dimensions  1 ’’Width  X4”  Height.  Additional  sound  absorbing  foam 
was  placed  around  the  sample  to  ensure  that  no  sound  was  leaking  around  the  sample 
particularly  at  steeper  angles.  Transmission  measurements  through  the  grating  were 
normalized  with  measurements  without  the  sample  in  place.  Each  measurement  was 
repeated  several  times  to  ensure  repeatability  and  consistency. 


3.  Numerical  Calculations.  The  numerical  calculation  was  carried  out  by  finite 
element  methods,  as  implemented  in  the  COMSOL  Multiphysics  software  package. 
We  assumed  the  sound  wave  impinging  on  the  sample  as  incidence  plane  waves. 
Periodic  boundary  conditions  were  set  around  a  unit  cell,  while  perfectly  matched 
layer  (PML)  absorbing  boundary  conditions  were  used  at  the  top  and  bottom 
boundaries  of  the  cell.  The  transmission  was  computed  by  a  ratio  of  an 
integration  of  input  and  output  pressure  field  at  the  top  and  bottom  boundaries  before 
PMLs. 
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